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The same principle is our guide in the placing of the products 
here. For the first product we have 3 + 4=7, and for the second 
3 + 6=9, the latter falling two places in advance of the former, in 
consequence of the occurrence of a cipher. For the next we have 
12 — 7 = 5; that is, the fifth decimal place of the second multi- 
plicand is the last effective figure, its product falling in the twelfth 
place. And here, also, the operation changes into that of con- 
tracted multiplication. 

It is with a view to perspicuity that I have exemplified the 
three processes, as each of the last two follows from that which 
precedes it by a natural sequence. But it will be understood that 
it is the tables for the application of the three-figure process alone 
that accompany these papers. The superiority of the process last 
named over the others, in regard to the number of tabular entries 
and of figures employed, is equally apparent in the anti-logarithmic 
as in the logarithmic application. It will be still more fully 
appreciated after the trial of a few examples. 

I must postpone till another opportunity a description of the 
methods employed for the construction and verification of the 
tables, and also an account of the origin of the processes that 
have been exemplified. 



On Interpolation, Summation, and the Adjustment of Nuinerical 
Tables. {Part III.) By W. S. B. Woolhouse, F.R.A.S., 
F.S.S., Vice-President of the Institute of Actuaries, Sfc. 
[Read before the Institute, 27th February, 1865.] 
The Adjustment of Numerical Tables. 

W HEN a series of quantities which depend on a fixed law, whether 
known or implied, follow each other in a due order of succession, 
the general accuracy of their numerical values may be satisfactorily 
tested by observing the regularity of the progression of a suitable 
order of differences. If the tabular quantities are the results of 
calculation from a given formula, with equidistant arguments, by 
differencing them up to a certain order the existence of an isolated 
error, if one should exist, is thus prominently exhibited and there- 
fore speedily detected. Moreover the precise locality of the error 
and the fact of its individuality being indicated by a characteristic 
interruption of the law of progression, the proper correction is 
ascertained by carefully noting the central position of the dis- 
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turbance, and revising the corresponding portion of the original 
computation. This work of revision may, however, be dispensed 
with, if necessary or desirable, as the magnitude as well as the 
locality of the indicated error can in general be readily deduced, 
with sufficient accuracy, directly from the differences, by the 
method of adjustment here given. The method, being so far 
independent of fundamental origin, is therefore equally applicable 
to the adjustment of a series of quantities derived from observation, 
or from statistical data ; only in this latter case, as almost every 
quantity must be more or less affected by imperfections arising 
from defective information and other known or unknown incidental 
causes, the influence of which may sometimes be augmented by 
paucity of numbers, the progression of the differences is disturbed, 
by a more complicated combination of the effects of au irregular 
series of errors, and it becomes more difficult to eliminate the 
value of each separate correction. 

In constructing and adjusting numerical tables, more especially 
those drawn from statistics, the method commonly employed is 
simply that of taking as data only a few of the unadjusted tabular 
values, the same being distributed apart at certain assigned con- 
venient intervals, and thence determining all the intermediate 
values by a suitable process of interpolation, so as to impart to the 
results a desirable graduated progression. Although expedient as 
regards facility of calculation, this practice is in other respects 
devoid of scientific principle. The meagre data upon which every- 
thing is made to depend may indeed consist of values amongst 
those most in error, and under any circumstances an undue and 
exclusive presumption of accuracy is tacitly given to a few values 
arbitrarily chosen to the exclusion of all the others ; and by using 
these values without any preliminary adjustment, the entire table is 
necessarily and systematically vitiated by whatever imperfections 
they may happen to possess. The real problem to be solved is to 
determine a method of constructing a properly graduated table, 
which shall represent the whole of a given series of values, or a 
given mass of original facts, with a greater degree of accuracy than 
that of any other table of a similar kind. A thorough investiga- 
tion of the subject, thus considered, and taken in all its bearings, 
may be said to open out an entirely new branch of science ; and 
in the discussion given in the present paper it will, I think, be 
admitted that some progress has been made towards its mathe- 
matical development. 

The uniform interval between the consecutive arguments of the 
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table being supposed to be sufficiently small, it will follow that if 
the exact functions were differenced, the differences would rapidly 
converge and exhibit no sign of discrepancy. But the given 
functions or quantities involve, in addition to the exact functions, 
their respective errors; and it is evident that the amounts of 
variation of the several differences would be produced by separately 
differencing the series of errors which, on account of their irregu- 
larity and consequent divergence of the differences, soon become 
prominently conspicuous. 

Isolated Errors. 

First consider the case of a single error. Let the value V be 
in error so as to require a correction v, and let the values which 
immediately precede and follow it be supposed to be exact. 
Then the effects produced on the several orders of differences by 
the correction v will be as follows : — 
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These are simply deduced by successive differencing in an ex- 
panded form, the absence of a symbol and the insertion of hyphens 
being always accounted as zero ; and it will be observed that the 
numerical coefficients which appertain to each order n are those of 
the expansion of the binomial (1 — 1)". 

Retaining the notation before employed ("Interpolation," 
vol. xi., page 64), the corresponding portion of the differences of 
the given series of quantities will stand thus : — 
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Given 

Values. 
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Now the quantity V , when corrected, will be V + i> ; and the 
adjustments produced by v on these differences will in like manner 
be obtained by incorporating the respective values contained in the 
former table. 

To detect the existence, as well as to determine the value of the 
correction v, it will be most convenient to have recourse to an even 
order of differences, as the greatest disturbance is then exhibited in 
the central difference, or that which stands opposite to V , the 
quantity in error. The practical methods of calculation which 
follow will now be obvious. 



Determination of an Error from Second Differences. 

An error in the value of a quantity V will be indicated by a 
marked disturbance of three second differences, A_j, b , b +1 , and 
the required correction is to be inferred so as to reduce them to an 
arithmetical progression. Thus, when the corrected second dif- 
ferences 5_ t +», b —2v, b +1 + v, are in arithmetical progression, the 
two third differences, C-i— 3v, c +l + 3v, derived from them, will be 
equal, and consequently the fourth difference, d +6v, will vanish 
or become neutralized ; 

.-. the correction («,) £ = ( b °- b -J + ( b °~J±}) . 

Hence the following simple rule, which may in general be applied 
mentally or by mere inspection, it being, of course, always under- 
stood that proper attention is given to the algebraic signs + and — . 

l 2 
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Rule. — From the middle second difference subtract separately 
the preceding and following differences; add together the two 
remainders, and divide by 6. The result, with its proper algebraic 
sign, will be the correction to be made in the value which stands 
opposite to the middle difference. 

Example. — To test the accuracy of a consecutive series of 
tabular numbers, they are differenced as follows : — 



Tabular 

Numbers. 


First 
Differences. 


Second 
Differences. 


Corrections of 

Second 

Differences. 


Adjusted 

Second 

Differences. 


1728 
2197 
2744 
3375 
•4109 
4913 
5832 
6859 
8000 
&c. 


+ 469 
547 
631 
734 
804 
919 
1027 
+ 1141 
&c. 


+ 78 
84 
103* 
70* 
115* 
108 
+ 114 
&c 


-13... t> 
+26... -2» 
-13... v 


+ 90 

96 

+ 102 



Here the three second differences distinguished by an asterisk are 
decidedly irregular as regards progression, and show that the 
tabular number opposite to the central difference ( + 70) must be m 
error. And for the adjustment of this error, according to the rule, 
we have 

middle— preceding diff. 
middle— following diff. 



Irregular 

Second 

Differences 


' +103 „„ 

70 ' 
— 4f> 
+ 115 40 * 


6J-78 


Oorrec 
V 


tion t' . . . . — 13 
4109 



Adjusted value 4096 
By replacing 4109 by this adjusted value, and revising the diffe- 
rences, they will now observe a regular progression, showing that 
no other error exists amongst the numbers put down. The second 
differences may, however, be otherwise directly and independently 
adjusted by applying to their original values the respective correc- 
tions + v, — 2i>, +v, which is already done at the right hand of 
the above table. 

The tabular numbers quoted in this example are the cubes of 

the numbers 12, 13, 14, 20. By differencing them out to 

fourth differences, even a single unit in error is at once disclosed ; 
thus : — 
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Tmre Valbks. 


Numbers. 


Differences. 


1728 
2197 
2744 
3375 
4096 
4913 
5832 
6859 
8000 


631 on 6 o 
721 Z 6 " 

817 102 6 
919 na 6 n 

+ 1141 + 114 



Primary Ebsob or an Unit. 


Numbers. 


Differences. 


1728 
2197 
2744 
3375 
•4097 
4913 
5832 
6859 
8000 


+ 469 j. -7R 

547+ r? + 6 , , 
631 \\ 7 + 1 
722 qi 3,* 
816 £* 9 + « 

1027 }5« + 6 + 1 
+ 1141 + 11 * 



The differences ( + 1, —4, +6, —4, +1) here agree in cha- 
racter with the fourth difference column of the first table, and 
thereby show that the indicated error is the only one that exists. 

Determination op an Error from Fourth Differences. 

The value of the correction of a quantity may in like manner 
be deduced from three disturbed fourth differences, so as to fulfil 
the condition that the three differences when adjusted shall form 
an arithmetical progression. The corrected fourth differences are 
here </_,— 4», <^+6», rf +1 — 4»; and when these are in arithmetical 
progression, the two fifth differences, e_ l + 10w, e +i — 10p, derived 
from them, will be equal, and therefore the sixth difference, 
/ — 20t>, will vanish or become neutralized; 

,. the corrects (.)= ^ = (*■-*>+(**-*> , 

which suggests the following rule : — 

Rule. — From the preceding and following fourth differences 
separately subtract the middle difference; add together the two 
remainders, and divide by 20. The quotient, with its proper alge- 
braic sign, will be the required correction, and will be applicable to 
the quantity which stands opposite to the middle difference. 

By the foregoing table, showing the effects produced by v on 
the differences, it will appear that five fourth differences will be 
subjected to alteration, and that they may be adjusted, with respect 
to v, without redifferencing, by directly applying to them the 
several corrections +v, — 4t>, + 6p, — 4i>, +». If the differences 
so corrected, when taken in connexion with those that immediately 
precede and follow them, should then be found to observe a satis- 
factory progression, the fact will show that the primitive values 
immediately adjoining V are free from material error. But if any 
marked irregularities should yet remain, other corrections will 
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require to be made in a similar manner ; and by thus continuing 
to repeat the process at points of greatest derangement, all the more 
prominent disturbing inequalities of the given series of numbers 
may be ultimately regulated. 

Whether second or fourth differences should be used in these 
calculations for the adjustment of any proposed set of numbers, 
must, to a certain extent, be left to the practical judgment of the 
computer. Perhaps fourth differences will be found to be the 
more generally applicable. As a general guide, however, it may 
be advisable first to consider what even order of differences of 
exact functions ought to be small enough to be practically disre- 
garded, and then to employ differences of the next less even order. 
And here it will not pass unnoticed that the foregoing simple 
rules of calculation may, if required, be readily adapted to higher 
orders of differences. 

The rules, as they are here given, were, many years ago, uni- 
formly and successfully employed by myself in officially examining 
and testing the accuracy of the results of the various and extensive 
computations made in the Nautical Almanac Office; and I also 
availed myself of their application in the adjustment of both the 
"Indian Military" and the "Experience" tables of mortality.* 
In the construction of all kinds of tables these rules will in general 
be found to be very effective in the adjustment of isolated or pro- 
minent irregularities, as they proceed upon the sound principle of 
applying each successive correction at the very point where the 
probable or certain existence of an error is most strongly indi- 
cated. If, for example, a number of errors, each of which may be 
of any magnitude, were promiscuously introduced into a page of a 
table of logarithms, or in any table computed mathematically 
according to a given formula or generic law which is still inherent 
to the general mass of values, the preceding rules, when applied to 
the differences, would necessarily have the effect of practically 
eliminating the several errors and restoring the table to its original 
state of accuracy. In all ordinary calculations involving the con- 
struction and adjustment of tables what has already been laid down 
may therefore be considered to meet every necessary requirement. 

* « Investigation of Mortality in the Indian Army. By W. S. B. Woolhouse, 
F.R.A.S., &c.," and " Tables Exhibiting the Law of Mortality deduced from $e Combined 
Experience of Seventeen Life Assurance Offices" It is presumed that no existing gradu- 
ated table of mortality represents so closely the general body of facts on which it is 
founded as the tables here referred to. Mr. Jellicoe, the President of the Institute, at 
the end of a valuable paper inserted in vol. i., page 166, has given a D and N table, 
together with the values of an annuity and annual premium, calculated from the adjusted 
Bengal Military Mortality, at 4 per cent. 
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Multiform Errors. 

We have yet, however, to enter upon a discussion of another 
part of the subject, which, if it be less essential in ordinary practice, 
nevertheless admits of most important applications, and also pos- 
sesses considerable interest in a mathematical point of view. If 
the generic law be unknown, and the table — as in the case of a 
table of mortality — be drawn from observation, then no part of it 
can be predicated as really or even probably accurate, and therefore, 
after the more prominent displacements and straggling errors have 
been adjusted as before described, it may afterwards be requisite, 
if a continuous and accurately-finished table be desired, to have 
recourse to some other more recondite process in order to finally 
eradicate a consecutive series of lesser imperfections, the disturbing 
effects of which on the differences are of course blended together 
without exhibiting any individual traces. 

To proceed to an investigation of the subject in this more 
complicated form, we have now to discuss the question of the com- 
plete adjustment of a combination of multiform small errors, or 
the neutralization of the assigned discrepancies of an irregular 
series of differences of a given order. 

In the first place we shall consider the case in relation to fourth 
differences. Let the given original quantities be denoted by 
V , V u V 2 , . . . . V„, and the sought corrections, with their proper 
algebraic signs, respectively by v , »„ v 2 , . . . . v n ; and suppose 
them to be separately differenced as follows : — 



V 
V 2 

v 4 
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7i 


«3 


A. 


72 


«4 








A.-2 


7n-3 


"»-l 







<$d 



R-4 



Then V -f » , Vj + r;,, . . . . ~V„ + v„ are the corrected functions, 
and their corresponding fourth differences will evidently be 
^o + S^O) d l + Bd l , .... d m .. i + Bd m _ i , where 8 is employed as a 
characteristic of variation, and not as a quantity. In most cases 
it may be practicable to assign such small values to Sd 0) Bd u 
.... &4-4 as shall impart a tolerable progression to these cor- 
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rected differences, and thus neutralize the principal irregulai-ities 
of the given differences d , d l} . . . . rf„_ 4 . The values $d , Bd t , 
&c, thus found, will be the only data for the determination of 
v , t>i, . . . . v n ; and as the number of these data, or mathematical 
conditions, is here less by four than the number of unknown 
values, the problem remains so far essentially indeterminate. This 
is indeed evident from the consideration that by arbitrarily assuming 
any four numbers for the unknown leading differences v , a , 
/3 , 70, and using these numbers in combination with the assigned 
fourth differences Sd , Sd u . . . . <5d„_ 4 , the several values of v, in 
common with the whole scheme of differences, become directly 
deducible by summation. The gist of the present inquiry is to 
assign such values to v , a , j3 , jo> tnat tne resulting displace- 
ments, v , »i, . . v„, of the original functions, shall be as small as 
possible. To effect this object the general amount of disturbance 
is most satisfactorily estimated by adding together the squares of 
the several corrections, or errors, for by that means positive and 
negative disturbances, which are alike objectionable, become all of 
them positive, and therefore accumulative. Hence the problem, 
mathematically enunciated, resolves itself into the following : — 

Problem. — To find a series of small quantities, v , v u . . . . v n , 
which shall have a given set of fourth differences, and such that the 
sum of their squares shall be a minimum. 

Before proceeding with the investigation of this problem, it 
may not be out of place here to premise a brief outline of the 
practical resolution of equations of condition by Gauss's " method 
of least squares," as the same may not be generally known or 
understood. At least, I believe it has not before been explained 
on elementary algebraical principles. 

Let a series of errors be expressed by the linear equations 
e =a a;+%+c z &c. 
e l =aix+b l y+c 1 z &c. 
e 2 =a2X+b^+C2Z &c. 
&c. &c. 

where « . &o > c o > & c - are given coefficients; let the number of 
equations exceed the number of variables x, y, z, &c. ; and let it 
be required to find values of these variables, such that the sum of 
the squares of the errors, or 

shall be a minimum. 

When this last condition subsists, any small alteration, whether 
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positive or negative, in the values of one or more of the variables, 
must increase the value of S. Suppose x to become x + e; then 
the series of errors will evidently become e + a e, ei + «is, e 2 + a 2 s, 
&c, and therefore the sum of their squares will be 

S+ 2(^0 +«,«,+ a^s, )e+(al+a\+al > 2 . 

Here the third term, involving c 2 , being essentially positive, it 
is evident that the value of the expression will always exceed S, 
whether s be positive or negative, provided that 

«0«0+«l<'l + a2«2 &c. = 0. 

And in a similar manner, by separately varying y, z, &c, it may 
be inferred that the same condition will be further maintained 
when 

Vo+Vl + &202 &C.=0 
C e + C l e l + C 2 e 2 &C.=0 

&c. &c. 

If in each of these equations the errors e , e t , e 2 , See., be 
replaced respectively by the preceding expressions, they will serve 
for the determination of the variables. Hence the following rule: — 

Rule. — Multiply the equations severally by the respective coeffi- 
cients of one particular variable, and equate the sum of the products 
with zero. By proceeding thus with respect to each separate variable, 
as many equations as variables are obtained ; and by solving these 
equations, of the first degree, the values of the several variables are 
ultimately determined. 

To apply this method of least squares to the present inquiry it 
will be requisite to express the general value (v m ) in terms of 
v o> «o» /3o> 7o» an ^ * ne given fourth differences Bd 0) Sdi, .... Sd„- t . 
This will be readily effected by summation. Thus, we have 
evidently 

7o =yo 

yi =y<>+&4 

Yi =y +H+M 

ya =y +Sd B +Bd l +Sd, 

Vm-s= , Yo+tid <) +8d l + $d i . . . . -\-Hd m _, 
.•.A»- S =A>+2y 

=/3„+(w— 2)y +(m— 3)H+(»»— 4)H &c. 
=/? +(TO-2)y +Z,{(w— « + 8)3<*,}. 
And, bearing in mind, the usual formula for the summation of a 
common factorial, viz., 
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* *?+!)■... (p+n)= KP + 1) - n+ i P+n+1) ..-(/) 

we similarly find 

n . , (to— 2)(m-l) 



^ (rn- x +3)(m-x+2) , 



and /. » m =» +So 



(to — 1W . (to — 2Yto— 1)to 
=t> +m«o + V 2 ' &,+ • f^ - 70 



2.3 

where 2 indicates a summation with respect to x from to m— 4. 

When m takes the consecutive values 0, 1, 2 ... n, the formula 

(8) will express with mathematical accuracy the several values of 

the series of corrections, and as the coefficients of the four variables 

„ x- i i m(m— 1) m{m— l)(m— 2) 
«\» «o» Po> 7o> a™ respectively 1, m, —^-g — '-, -+ ^ '-, 

the resolving equations, by the method of least squares, so as to 
make Sc 2 a minimum, will be, according to the foregoing rule, 

X _ ~ „♦»(»» — 1) * „m(m— l)(m — 2) 
2i> m =0, 2TOt> m =0, S v - ' v n =0, S-i ^ 'v m =0; 

and these are evidently equivalent to 

2» m =0, Sto»„=0, Sw a p M =0, Sm 3 « M =0 .... (e). 

In effecting these summations with respect to m, from to n, 
after substituting the formula (8), the term involving §d x will, in 
each case, consist of factorials, by further substituting the alge- 
braical equivalents, 

m =(m—x)+x 



m 2 =(m— x)(m— x— \ ) + {2x— 1).(to— x)+x* 

m s =(m— x)(m— x— l)(m— x— 2) + 3(x— 1) . (m— x)(m— x— 1) 
+(3« 2 — 3x + l).(m— x)+x*. 

To abbreviate the process, let 



(n—x+3)(n—x+2)(n—x+l)(n—x). 
— Id, . . . ( 9 ); 



also, 
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M=2 / , s 
P=2/i,{4»+0M-4)} 

R=2^{20» 8 +(10tf+50)n a +(4* 2 + 26:r+38)» 

+ (s 3 +8s !! +18:H-8)}. 

Then, for all consecutive values of m from to n, the results of the 
respective summations contained in (e) will be 

/., 1V (»+!> ■ (n + l)n(n-l) 
(»+1>oH g — ^ 23 ^° 

(n+lXn-lX — 2) „ . M _ 
+ ~%U ro+ 4 ~°' 

(»+l)n . (» + l)n(2n+l) (n + l)«(n-lX3»+2) . 
— 2 — »o+ ^ «o + g^ ft 

( w +l>0»-l)(n-2X4n+3) P 

+ — 2JZ~5 ro+ 4^ -°' 

(»+l)n(2»+l) (»+l) 2 » 2 , (»+l)n(«-l)(12» 2 + 15»+2). 
£3 »o+ — «o + 2^X5 A> 

(»+l)»(»-lX»-2X10»* + 14n+3) , Q n 
+ 3X5."6 ro+ 3X5 -°' 

(n+l)*» s (»+l)«(2n+lX8« 2 + 8»— 1) 
__ , 0+ _ « 

(n+l><w— l)(10n s +18n 2 +5n— 2) 
+ 2^X5 Po 

(n+l)»(»-lXn-2X20n 8 +40n 11 +16»-3) R _ n 

+ 4.5^7 yo+ 4X?- ' 

Dividing these equations by the respective coefficients of v , they 
become 

n , nin-l). , n(n-l)(n-2) M 

"° + 2 ° 0+ ~2T~ <* 0+ 2T3l ro+ 40hT) = °' 

l 2»+1 . (»-lX8n+2) fl 

c ° — 3" — "° Ha 

, (»-lX»-2X4»+3) , , P_ 

+ "3X5 yo+ 10n(»+l) ' 

3 n(n+l ) (w-lX12n 2 +15w+ 2) . 

" 0+ 2(2.+ 1) a ° + 43(2^+1) ^ 

(n-lX«-2)(10w 2 + 14»+3) Q_ _ Q 

+ 3.4.5(2n+l) 7o+ 10»(n+lX2n + l) ' 
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2(2n+l)(3ra 2 +3n-l) (rc-l)(10n 3 +18n 2 + 5«-2) n 

V ° + 15»(«+1) Oo+ 30«(«+l) ft 

(n— l)(n— 2)(20» 8 +40» ii + 16»— 3) R_ 

+ 210»(»+1) r ° + 35w 2 (n+l) 2 

Differencing these equations, and then dividing by the respective 
coefficients of a Q , we get 

, n-1 „ , 3(n-l)(n— 2) 15rcM-6P 

«o+ -5- M ™ Yo=~ 



2 r " 20 '" 10ra(»+l)(»+2) 

, 3»-l fl , »(»-2) 6(2n+l)P-6Q 

a ° + T~ A+ ~5~ yo_ 10(«-l)....(n+2) 

4» 8 + 3» 2 -w- 2 r (»— 2)(10» 8 +12» 2 +w— 3) 
a ° + .2(3n 2 +3»+2y^ 0+ " l4(3» a + 8»+2) 7o 

21w(w+l)Q-6(2w+l)R 



,7(»— 1). . . .(«+2)(3« 2 +3w+2) 
Again differencing, and dividing by the coefficients of /3 , 
„ »-2 15M 18»P— 6Q 

Po+ -5-yo=» 



2 '° (m+1)(»+2)(»+3) ' (n-l)....(n+3) 

60R 



a 4n-5 __ 6(3« !! + 3»+2)P-12(2w+l)Q 
Po+ ? 7o- (n-2)....(«+3) 



7o= 



7(»-2)...(»+3) 
210M 84(6n 2 -3n+2)P-420«Q + 120R 



(n+l)....(n+4) (n-2)....(n+4) 



30(4w-5)M 6(45n 2 -9ra+ 14)P— 24(9» + 1)Q+ 60R 

^o - (n+l)....(»+4) + (»-l)....(«+4) 

30( W JJ_ W+ 1)M 6(20w 2 + 5w+8)P-6(15»+4)Q+24R 

a ° — (n+1) («+4) ~~ 2.» (»+4) 

_ 2(2w+l)(» 2 +n+3)M— (6w 2 + 6»+5)P+2(2n+l)Q— R 
* 0_ (n+l)....(»+4) 

Hence, by restoring the foregoing values of M, P, Q, R, we find 

_ ~ (»+lX«+*X«+8) 
V °~ i ( w+ l)(„-|-2)(n+3)(n+4)'' 1 

(*+l)(*+2){n+i(x+4)} 

a °~ «(»+ 1 X»+ 2 X»+3)("+ 4 ) 



, (x+ l){»g +(4tt+15>+(10:r 2 + 66:r+104)} 
(n— 1)» (n+4) 



Po= 6S 1, i\„ 7TTT\ /*» 



__ tt3 + (4s+13)n 2 +(10^+58a;+74>H-(20^+16()s !i 4-388*+272) 
yo- ~ 6S («-2X»-l>. ■ • .(»+4) ^ 

These results are perhaps more systematic when put in the 
following form : — 
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p(*-H)(s-H2)0+3) 



*° S (»+!)....(«+ 4) 



a °-- 3S (»+l)....( W +4)(, 1+4 — J 

A.- 6S ohTi). . . . (w+4 )l, 1 + 4 ir +10 -k,^i)-J 

7o= _ 6s __JL__f 1+4 =±* + i C?+*Xf+8) 



(n+l)....(«+4)V » n(»»— 1) 

+2 ° «(»-l)(«-2T J ' ' " * ( } ' 

From the values of »q> a , /3 , y , thus determined, the several 
values of v can be successively deduced by summation according to 
the scheme of differences, or they may be calculated separately and 
independently by substitution in the formula (2). 

To thus extend the investigation to the adjustment of higher 
orders of differences would involve the discussion of equations which 
increase both in number and complexity, and would evidently be a 
work of considerable labour. I have, however, by a somewhat 
novel artifice, succeeded in reducing the fundamental expressions 
to a form so remarkably unique and comprehensive, that the 
general solution of the problem, for any order of differences, 
becomes not merely simplified but divested of all complication 
whatever. 

The primary basis of this improvement consists in prefixing 
to the scheme a certain set of supplementary differences obtained 
by a receding process in which the prefixed values of the 
original function are severally assumed to be zero. These sup- 
plementary differences will be placed in the form of a wedge, 
and will cause the differences of each order to be equal in number 
to that of the original functions, and impart to them the pecu- 
liar property of making the values of the several functions depend 
exclusively on any stated order of differences according to a formula 
in which the coefficients are all included in one common factorial 
expression. 

Thus, resuming the case of fourth differences, and retaining 
an uniformity of notation, the extended scheme will become as 
follows, where the upper ridge of values v , a- 1} /3_ 2 , y- 3 , &?_ 4> 
are necessarily equal, though designated by different symbols : — 
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*>3 
«6 



<*a 

«3 
«4 



/3-a 


y-8 


P-i 


y-a 


A> 


y-i 


A 


yo 


A 


yi 


A 


ya 


• 


y»-s 


&-« 





oa_ 3 



W 



Also the formula (8), when made to commence with S«L 4 , so 
as to include the four supplemental differences, will become simply 



(m— *+3X»»— «+2)(»»— ar+1) . 



.(e). 



For the coefficient of Bd x being a function of m— x, is independent 
of the epoch ; and when the epoch is transferred to »_ 4 , the values 
of v , a , /3 , y , are 0_ 4 =O, a_ 4 =0, /3_ 4 =0, y_ 4 =0. The 
terms comprehended under 2' (accentuated) will now begin with 
x— —4 and end with x=m— 4. 

When SpJL is a minimum, it has already been shown, equa- 
tions (e), that 

2» M =0, 2mv M =0, 2m 2 e m =0, Sm 8 f m =0. 
And as these, when compounded with any finite functions of x, 
must still vanish, they are evidently equivalent to the following : — 

2v m =0 

S»„(»j— *)=0 __ 

2v m (m—x)(m— x— 1 ) =0 

Sv m (m— stym— x—l)(m— x— 2)=0. 

When the value of v m , by (*), is substituted, the several 
coefficients of $d x are now of the common factorial form (/) ; and 
if, as before, we abbreviate by putting 



(n— x+3)(n — x+2)(n— x+l)(n — x) , 



is) 



the summations, after being taken with respect to all consecutive 
values of m from to n, will give 

2>,=0 

Zpj(n— x— 1 )=0 

S>,(»— x— !)(«— ar-2)=0 



2>,(w— «—!)(»— *— 2)(n— «— 3)=0. 
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And these are again equivalent to 

2>« =0 < 

2 V* =° (6). 

2V/i,=0 , 

Hence, when 2 (unaccentuated) does not include values which 
depend on the four supplementary differences, these last equations 
will give 

2/t„ =— /*_ 4 — p-3— P-2— /*-i 
Ssju, = 4/j_ 4 + S/*_3+2/i_ 2 +/i_i 
2*V»=— 16/i_ 4 — 9jk-3— 4ju_a— /*_i 
2afy,= 64/i_ 4 + 27/<_3+8/i_ 2 +/i_ 1 

for the determination of the four unknown values ju_ 4 , /u_ 3 , /t_ a , 
ju_i, and thence 8<?_ 4 , $d- 3) &L 2 , 8rf_j. 

Adding each respective equation to the next following one, 

we get 

2(*+l)/x, = 3/i_ 4 + 2/i_3+ /*_ 2 

S(*+l>/«, =-12jU_4- 6/l_3-2/l_ 2 

J5(*+l>fy,= 48/t_ 4 +18/i„3+4/t- 2 ; 
and, adding twice each of these to the next following, 

2(s+l)(*+2K =-6/i_4-2/*_ 3 
2(ar+l)(ar+2>/i,= 24/i_ 4 +6/«_ 3 ; 

•' i - 4 = S 2^ * 

^ (s+l)QH-2)(*+4) 
/*.»=— 2 5 {!*■ 



Similarly we find 

/*-2 



(x+lX f +3X f + 4) 



„ (»+gX«+3X«+*) 
/«-i=- S 5"3 



2.3 *" 



But, as assumed in (g), 



( W +l)(n+2)( W + 3)(n+4) SJ 
^- 4= 23 M ~* 

n(»+l)(» + 2)(»+3) s , 
/*-s= g^ M - 3 

(»-lK»+l)(»+2) s , 
/i. 2 = ^ <5rf- a 

( w _2)(»-l)n(»+l) SJ 
^ 1= 2^ **-> ' 
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(«+l)(«+2X*+8 ) 
(» + l)(„ + 2)(»+3X«+4)' 
(«+lX«+2X *+4) 
n(n+lX»+2X»+3)' < * 

. S («+1X*+»X*+*) 
(»-l)«(n + lX»+2) / '* 

__ (^+2X^+3X^+4) 

(n_2X»— !)»(«+ 1)''* 



. (A), 



in which the summations under 2 are from x=0 to *==»— 4 
inclusive, and the complete form is embodied in each term. 

By algebraically differencing back from the symbols v , a , 
fioi 7o» the algebraic values of the supplementary differences are 
found as follows : — 



n 



»0 
"0 



»0 



— 2» +oo 
»o— «o+0o 

ro 

ao=3M.4 + W_3 
j3o=88rf_4+2Jrf. 8 + 3rf_ 2 



#0 =3«?_4 

— 3f +«0 =&L 8 

3t> — 2a +/3 =8rf_ 9 

— «-'o+«o— Ai+yo =8rf_i 



(0, 



which, after substituting the values (k), will be found to be equiva- 
lent to the formulae (h), before determined. 

The generalization of the foregoing method, for any order (p) 
of differences, is perfectly analogous to what has been given, and 
the systematic law which obtains in the form of the final results is 
such that they may be readily put down for any other order 
of differences. It will be observed, that in the expressions for 
Bd-p, . . . §d_i, the numerical coefficients are the respective terms 
of the expansion of the binomial (1 — l)*" 1 ; and that in each 
numerator the absent factor corresponds with the subnumeral 
indice of Sd. 

Thus, for sixth differences, we should have 

w (ot— x+5)(m— ~x~+ 4)(m— x + 3)(m— x+2)(m—x+i) s , 
r »= S 2XL5 J " "^ >' 

(n— i+5X»— «+4X w_ «+8X»— «+2X»-«+lX»»- a: )s f , ... 
fV= 2^X5 d - / " ^ )5 



and, for the supplemental differences, we should obtain 
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(»+lXa>+2)(g+3X*+*X«+S) 

(«+l)(«+2)(n + 3)(« + 4)(«+5)(n+6) / '* 

(*+lX«+2X*+8X*+4X»+6) 
n(n+ l)O+2)(n + 3)(»+4)0i + 5) ^ 

102 ^+1)^+2)^+3)^+5)^+6) 
(n-l)n(»+l)(»+2X»+3X»+4) P * 

_ 102 (*+lX*+2X*+4X*+5)(*+6) „ 
(n-2Xn-l)»(n+i)(»+2)(» + 3)^* 

5Z Qe+1)0+3)(*+4)Qe+5X*+6) 

( w _3Xn-2)(n-l)w(w+l)(»+2) /ia: 

(^ + 2)^+3)^ + 4X^+5X^+6) 

(n-4)(n-3X»—2X»— !)«(» + 1)'*'* ' ' K } ' 



Practical Adaptation of Fourth Difference Formulae. 
To simplify the formulae for calculation, let x'=(n— 4) — x, and 
tabulate the factorial functions 

(*+l)(*'+2)(*+3)0r'+4) 



k= 



2.3.4 



, 4*. 12*. , 12*, 



x+l' 



Then we shall have 



t kp 



*- 3 =2 



kg 



also, 
that is, 



ScL 2 =2— —ld x 

"n-2 

ks 
arf_,=2- — Stf, 

*n-3 



(»); 



«i= 4Sc?_ 4 — f- a</_ 3 

» 2 =10Srf_4+ 43rf_3+ Srf_ 2 

t 3 =20M_ 4 +10W_3+ 45rf. 2 + Jrf_, 
«4=858rf_4+208rf. 8 +108rf_j+4arf_ I +arf t) 

&c. &c. 

But it will be found that the several values of v may be otherwise 
readily deduced by successive summation of the columns of 
differences. 
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The tables here required are easily prepared, and are as 
follows : — 

Subsidiary Tables for Fourth Differences. 



X 


k 




X 


P 


? 


r 


S 


Sec. 


&c. 





+ 1 


- 4 


+ 6 


- 4 


20 


10626 




1 


4 


15 


20 


10 


19 


88SS 




2 


10 


36 


45 


20 


18 


7315 




3 


20 


70 


84 


35 


17 


5985 




4 


35 


120 


140 


56 


16 


4845 




5 


56 


189 


216 


84 


15 


3876 




6 


84 


280 


315 


120 


14 


3060 




7 


120 


396 


440 


165 


13 


2380 




8 


165 


540 


594 


220 


12 


1820 




9 


220 


715 


780 


286 


11 


1365 




10 


286 


924 


1001 


364 


10 


1001 




11 


364 


1170 


1260 


455 


9 


715 




12 


455 


1456 


1560 


560 


8 


495 




13 


560 


1785 


1904 


680 


7 


330 




14 


680 


2160 


2295 


816 


6 


210 




15 


816 


2584 


2736 


969 


S 


126 




16 


969 


3060 


2823 


1140 


4 


70 




&c. 


&c. 


&c. 


&c. 


&c. 


3 

2 


35 

15 


























1 


5 















1 















These tables may be applied with the greatest facility to the 
resolution of any given case by means of the following simple rule, 
suggested by the formula? (m) : — 

Rule. — Beginning with af=n, transcribe the several numbers 
to the end of Column (k). The first four of these will be the four 
divisors of the respective formula, and the remaining numbers will be 
the multipliers of each set. Opposite to these multipliers write down 
as many leading values of p, q, r, s, in four adjoining columns. 
Then the numerical coefficients of the formula expressing the four 
supplementary fourth differences will be found by taking the horizontal 
products kp, kg, kr, ks, and dividing the sets by the respective divisors 
before mentioned. 

Example 1. — Let there be only one fourth difference ; then 
w=4, and there will be only one value of x, viz., x=x'=0. 

The numbers from the tables will stand thus : — 



k\ 



I 70 
35 
15 


...k 


n 






5 


P 


9 


r 


s 


, 1 


+ 1 


-4 


+ 6 


—4 
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•'• W_ 4 =+ — W 

4 

W_ 3 =— — W 



6 2 

W_ 2 = + tt W = + T W 
15 5 



w_, 



5 H - 



Hence, multiplying by 70 to obviate fractions, tbe values of v are 
found by continued summation of each preceding column, as fol- 
lows ; the supplementary differences being cut off by horizontal 
lines : — 



70W 


70y 


70/3 


70a 


70w 


+ w 


+ W 


+ Wo 


+ H 


+ W 


— 8W 


- 7W 


— GBd 


— bM 


-4H 


+28W 


+21W 


+ 15W 


+ 10W 


+ 6W 


— 56W 


-35W 


-20W 


-1024 


-42rf 


+ 70Sd 


+35W 


+lbhd„ 


+ bSd 


+ W 



.-. r =+— W 



«'2=+ to W 
" 3= ~"70^ 



70 



Example 2. — Let there be two fourth differences, d , rfj ; then 
n—5, and the tabular numbers are 



126 










70 










35 










15 


P 


9 


r 


* 


5 


+ 1 


- 4 


+ 6 


— 4 


1 


4 


15 


20 


10 



.-. W_ 4 = 



126 H+ 126' 



4 3 

W_ 3 = — — W> _ TT Wi 



Sd_ 9 == 



-W + -Sd x 



4 2 

W_,= — - 24— g W, 



And, operating only with the numerical coefficients, 



1262^ 


126y 


126/3 


126a 


126« 


+ 5+4 


+ 5+4 


+ 5+4 


+ 5+4 


+ 5+4 


- 36-27 


- 31-23 


-26-19 


-21 — 15 


— 16 — 11 


+ 108 + 72 


+ 77 + 49 


+51+30 


+30+15 


+ 14+ 4 


—168-84 


— 91—35 


-40- 5 


-10+10 


+ 4 + 14 


126 


+ 35—35 


- 5-40 


-15-30 


-11-16 


126 


+ 35+91 


+30+51 


+ 15 + 21 


+ 4+5 



M 2 
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16 11 , 

*» 126 H " 126' 

14 ., 4 , 



11 jw 16 *,/ 



and 



Example 3. — Let there be three fourth differences ; then n=6, 



M 



f 210f 
126 










70 










35 


P 


9 


r 


* 


15 


+ 1 


— 4 


+ 6 


- 4 


5 


4 


15 


20 


10 


*• 1 


10 


36 


45 


20 



" 4 210 42 or 42 '^42 

„ kg 20., 25., 12, 

^= S I&=-42^-42^-42 < 



70 



14 



14 



14 



" l= "85" =_ Y °"~ T l ~ 7 2 
And, by summation of the coefficients, as before, 



42S<* 


42y 


42/3 


42a 


42w 


+ 3+4+2 
-20-25—12 
+ 54 + 60 + 27 
—72—60-24 


+ 3+4+2 
-17-21-10 
+ 37 + 39 + 17 


+ 3+4+2 
-14—17- 8 


+ 3+4+2 


+3 +4+2 
_8- 9-4 
+ 4 * -1 
+ 4+10+4 
-1 * +4 
_4_ 9-8 
+ 2+ 4 + 3 


-11-13- 6 

+ 12+ 9+ 3 

» +10+ 5 

- 5-10 * 

- 3- 9-12 
+ 6 + 13+11 


+ 23 + 22+ 9 
-12+ 1+ 2 
_ 5_20- 5 
+ 2+ 1-12 
+ 9 + 22 + 23 


_35_21- 7 
+ 7-21- 7 
+ 7 + 21- 7 
+ 7 + 21 + 35 


42 

42 

42 



v = -R+ -H+ -ia\ 



+ The numbers k put on a separate slip may be mechanically placed in position upon 
the other table, and thus all transcription will be dispensed with. 
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f 2 = 


5* * -5** 


l'3= 


i^o+j2 Wl+ 42 H 


»4=- 


"42^° * + 42 H 


« , s=- 


_ A H- & Wl -i2 H 



The accuracy of the work will in all cases be checked by the 
symmetry of the coefficients, as well as by the condition S»=0. 
Similarly for four fourth differences, or n=7, the results are 

*«*-«= ^K+^H+^4+ eg^S 

,, 28.. 45,, 36 14 

M -*= ~ 42 M °~ 42 Ml ~ 42 3 *~ 42 W »- 

2== 42 ° + 12 l + 42 2 42 3 

, 28 30 20 7 

&/_!=— — &? — 14 ' — 14 2_ 14 3 

and, by summation of the coefficients as before, 



f = 


7 12 10 4 


t',=- 


16,, 159,, 116,, 6 
"66 H -7T6 W, -7T6 W "~66 H 


t> 2 = 


4 40 59 4 


v 3 = 


^■S-^ + Ss^+S 8 * 


v 4 — 


|h+^- 6 ^ + 1^ + 4^ 


v & =- 


4 59 40 4 
" 66^-7X6^" 7^66 H+ 66 H 


«'6=- 


6 . , 116,, 159,, 16,. 
-66^-766^-7X6^-66^ 


v 7 = 


4^+^-M+§H+^^. 



Again, for live fourth differences, 
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.. 14 28., 30 s , 20 S7 , 7 

., 28., 35 18 35., 4 

10 45 2 

^.2= 22rf + — &*,+ — Srf 2 + 2M 3 +-M 4 

. 6 14 o 

20 25 50 25 4 
Id. , = — — Sd — — Bd 1 — — dd 2 - — Sd z - - Sdi ; 

14 28 30 20 7 s 
28 jh 91 ^ 42 ju 25 ^ 8 jw 

8 - , 22 . , 139 20 2 
t,6 =-99^-99^-l433^-99 M3+ 99 H 
8., 25., 42,, 91 28 

The formulae may in like manner be constructed for any pro- 
posed number of differences. And, finally, to adjust the given 
original values V , V,, V 2 , &c., the calculated corrections or 
adjustments v , »,, i> 2 , &c, are to be respectively applied according 
to their several algebraic signs. 

We have yet to make a further most important improvement 
and extension of the method, by the introduction of certain con- 
siderations, connected with the theory of the differences of alge- 
braic functions, which have the effect of concentrating the general 
principles of the investigation within the narrowest possible limits, 
and of finally bringing the subject under a remarkably comprehen- 
sive point of view. 

As the coefficients of ld x , in the equations (0), involve x alge- 
braically to powers not exceeding the seventh, it follows that the 
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eighth difference of any nine successive values of a set of each of 
these coefficients will be zero. Hence it will appear that these 
four equations will all of them be fulfilled if the several differences 
be replaced by any of the following horizontal sets of numerical 
values : — 



id_ 4 


id. 3 


id_ 2 


id , 


id 


idi 


id 2 


id 3 


id 4 


id, 


id 6 


id-j 




+ 1 


-8 


+ 28 


-56 


+ 70 


—56 


+ 28 


— 8 


+ 1 











&c. 





+ 1 


- 8 


+ 28 


-56 


+ 70 


-56 


+ 28 


- 8 


+ 1 








&c. 








+ 1 


- 8 


+ 28-56 


+ 70 


-56 


+ 28- 8 


+ 1 





&c. 











+ 1 


- 8[ + 28 


-56 


+ 70 


-56| + 28 


-8 


+1 


&c. 



&c. 



&c. 



&c. 



for each formula will then, in fact, express the eighth difference of 
nine consecutive coefficients ; and, for the reason just stated, each 
of these eighth differences of specified algebraic functions must 
necessarily vanish. Therefore, as the same must hold with any 
positive or negative multiple of a horizontal set of numbers, it will 
follow generally that all the requisite conditions (0) will be satisfied 
provided that 



°d-i= <po 

id. 3 =— 80 o + 0, 
id. 2 = + 280 o - 8^ + 
id_ l = —5S<)> + 2Sfi— 



P2+ 03 



id = + 70^0— 56^, + 28</> 2 — 80 3 
idi =—5600+700,-5602 + 2803 
id 2 = + 280o—560i + 7O0 2 -5603 
id 3 =— 800 + 280,-5602 + 7003 



• 0) 



(2) 



id 4 = 
id, = 



0o- 801 + 2802-5603 j 

• • • • 01 — 802 + 2803 
02 _ 803 



?3 I 



■ (3) 



where O > i>i> 02> i>3> represent a set of arbitrary or determinate 
quantities, which may evidently, in like manner, be extended to 
any number. These algebraic expressions are moreover the fully 
developed eighth differences of the series 

0, 0o» 0i) 02) 03» 0; 

and hence also the original quantities v , v x , v 2 , v 3 , v t , &c, are the 
fourth differences of the same series, and are therefore expressed as 
follows : — 
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V — 0o 

»i=-40 o + h 

»2=+ 600 40i+ 02 

»3=-40 O +60 1 — 40 2 + 03 
»4= 00—401 + 602—403 

t> 3 = .... 01 — 402+603 

^6= 02 — 403 

«7= 03 • • • • (4). 

The general equations, arising out of the principle of least 
squares, are thus effectually embodied in the simple condition that 
the required series of corrections (i>) shall be the expanded fourth 
differences of an interior and more concentrated series of values (0), 
four less in number, or the same in number as that of the fourth 
differences d , d u d 2 , &c. 

This view of the subject may be arrived at in a still more direct 
and simple manner by a consideration of the initial equations (e) 
in connexion with the fundamental formula (i). 

Suppose four inner sets of differences, denoted by w, y, z, <p, 
to be inserted in the scheme of adjustments, and to be determined 
from v , Vi, v 2 , v 3 , &c, by successive summation, thus: — 

M_ 4 



Foci 



<Po 



#- 











P-2 


dd_ 3 






»0 


a_i 


P-V 


Sd_z 




W 




"0 





y-> 5J 


y» 




Vl 




A> 


&Lj 


z 


K>1 




a i 




y ° TT 


y\ 




H 




Pi 


ca 


«i 


«>2 




a 2 




?> SJ 


y* 




*3 




Pi 


Sdi 


z 2 


M>3 




"3 


p 3 


72 • 


. Vs 




*>4 








«>4 




«4 




: WL-4 




2 


n 




A.-2 




*n~3 


od B _ 3 




W»-l 


»» 


«n- 


} P»-! 


2(4-2 








«» 


A. 


y " _l ML-. 
7» *, 



and by formulae analagous to (e) we shall have 
y n _,=Sw=S(n+ 1— *>. 



=i»-2=2y =2 



(n+1— s)(«+2— a) 



.=Es =2 



(n+1- 



2 
-*)(»+ 2- 
2.3 



-*)(n+8— «) 
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But, by equations (e), 2iV='2xv='2x i v='Ea?v=0 ; 

Hence these inner sets of differences will severally terminate 
with the values w„_i, y„- 2 , %*-%> ft>-4> as they appear in the 
scheme. If, however, any more inner sets of differences were to 
be constructed, it is evident that, unless S# 4 »=0, the number of 
terms, instead of continuing to contract, would become indefinitely 
prolonged. In consequence of this remarkable peculiarity, which 
is evidently applicable to any order of differences, the series of 
values of ^ are here designated the focal adjustments with respect 
to the differences under consideration. 

Whatever be the values of <p Q , <j> u <p 2) <j> 3 , &c, which form the 
radical basis of the differencing scheme, 2i> 2 will necessarily be the 
least possible out of all sets of quantities which can have for their 
fourth differences the series $d , Sd u $d 2 , &c, or the more extended 
series &L 4 , dd. 3 , §rf_ 2 , &c. Therefore, as this condition places no 
restriction upon the values of (j> , <j>i, <p 2 , <j> 3 , &c, they can be made 
to fulfil any other condition that may be required. The other sets 
of differences being linear functions of ft,, ft, <p 2 , <p s , &c, and 
exceeding them in number, are, however, not wholly independent, 
and it is evident that only as many conditions as focal corrections 
are admissible. 

When, for example, the requisite adjustments of four fourth 
differences are given, the corresponding focal adjustments are 
determined by the four equations (2), viz., 

70ft,— 56ft + 28ft— 8<t> 3 =$d 
— 56ft, + 70pi— 56ft+28ft,=H 

28ft,— 56ft, + 70ft— 56ft=&? 2 
— 8ft+28ft— 56ft, + 70ft,=&4 

which give 



(5); 



462ft =49H+ 84<Wi+ 70Srf 2 + 28Srf 3 \ 
462ft =84arf + 177Sd 1 + 164gtf 2 + 70W 3 
462ft = 70M + 1 64M, + 1 77H + 8iSd 3 
462ft,=28a<* + 70H+ 8ihd 2 +42Sd 3 



GO- 



The preceding and following supplementary differences are likewise 
found by substitution in (1) and (3), and the results will be found 
to accord with those of the same example on page 157. Hence also 
the focal values may be otherwise found by a continuation of the 
process of summation, with respect to the coefficients, thus : — 



162 



On Interpolation, Summation, §c. 



[April 



462b 


462i» 


462y 


49+ 84+ 70+ 28 


49+ 84+ 70 + 28 


49+ 84+ 70 + 28 


-112-159-116- 42 


-63- 75- 46-14 


-14+ 9+ 24 + 14 


28- 40- 59- 28 


-35-115-105-42 


-49-106- 81-28 


56 + 140+ 80+ 21 


21+ 25- 25-21 


-28- 81-106-49 


21+ 80+140+ 56 


42 + 105 + 115 + 35 


14+ 24+ 9-14 


-28- 59- 40+ 28 


14+ 46+ 75 + 63 


+ 28+ 70+ 84 + 49 


-42-116-159-112 


-28- 70- 84-49 


* 


28+ 70+ 84+ 49 


* 





462* 


4620 


49 + 84 + 70 + 28 


49+ 84+ 70 + 28 


35 + 93 + 94 + 42 


84 + 177 + 164 + 70 


-14-13+13 + 14 


70 + 164 + 177 + 84 


-42-94-93-35 


28+ 70+ 84 + 49 


-28-70-84-49 


* 


* 


Focal coefficients. 



Here the last set of numbers are the coefficients of the formulae ( p), 
and the symmetry of the numbers throughout is an effectual test 
of the accuracy of the work. 

The following series of tables, deduced from equations anala- 
gous to (5), exhibit the similar coefficients for any number of 
differences and foci up to seven. As they are all positive, the 
algebraic sign is not required. 

Tables A. — Coefficients for deducing Focal Values from their Eighth 
Differences or from assumed Fourth Difference Adjustments. 



□ 



126^ 






420 








4620 




5 

4 


4 
5 




3 


4 


2 




49 


84 


70 


28 


i 


7 


4 


84 

70 


177 


164 


70 






2 


4 


3 


164 


177 


84 












28 j 70 


84 


49 



1980 



42900 



28 


56 


60 


40 


14 


56 


133 


156 


110 


40 


60 


156 


204f 


156 


60 


40 


110 


156 


133 


56 


14 


40 


60 


56 


28 



756 


1680 


2100 


1800 


1050 


336 


1680 


4340 


5880 


5300 


3200 


1050 


2100 


5880 


8715 


8380 


5300 


1800 


1800 


5300 


8380 


8715 


5880 


2100 


1050 


3200 


5300 


5880 


4340 


1680 


336 


1050 


1800 


2100 


1680 


756 
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180 


432 


600 


600 


450 


240 


72 


432 


1188 


1776 I 1860 


1440 


786 


240 


600 


1776 


2868 


3176 


2560 


1440 


450 


600 


1860 


3176 


3743 


3176 


1860 


600 


450 


1440 


2560 


3176 


2868 


1776 


600 


240 


786 


1440 


1860 


1776 


1188 


432 


72 


240 


450 


600 


600 


432 


180 



The darker lines in each table show how symmetrically certain 
coefficients distribute themselves within the four sides of each con- 
centric rectangle. Also the several sets of coefficients are precisely 
the same whether they are read off vertically or horizontally. The 
symmetrical relations thus exhibited may indeed be inferred from a 
like symmetry in the coefficients of the original equations from which 
the tabular numbers are derived ; and they present an efficient and 
useful check on the accuracy of the calculation. An improved 
method of constructing these tables will be given hereafter. For 
more than seven values the resulting fractions are so cumbrous 
that it would be more convenient to tabulate them in decimals. 

If the necessary adjustments of a number of fourth differences 
can be first assigned, the focal adjustments may be calculated with 
these tabular coefficients according to formulae similar to (p), and 
thence the adjustments of the primitive quantities may be found 
either by successive differencing or by substitution in expressions 
of the same kind as (4). 

It will be observed that the foregoing methods of calculation 
proceed upon an assumed previous knowledge of the adjustments 
required to be made in the fourth differences, and that, in rectifying 
the progression, these must again depend in some measure upon 
the judgment and discretion of the computer. In order, however, 
to determine theoretically the values of $d , Sd t , $d 2 , &c, so as to 
effect the requisite adjustment of the quantities V with the least 
possible displacements, we may assume as an available new con- 
dition that the sum of the squares of all the amended fourth 
differences shall also be a minimum. Resuming the case of four 
foci, which will sufficiently indicate the general principles of the 
method which follows from this hypothesis, let the differences be 
carried out more extensively as follows : — 
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M>0 

y<> 

z w l 

fa yi 

Zl w 2 
fa Vi 

2 2 W 3 

fa 2/3 
Z 3 W 4 

fa y\ 

ys 



A-2 

u _l 

— £-1 
A 

A 
A 
A 
A 
A 



"7 



A 



y-3 
y-2 

Y-i 

yo 
yi 

72 
73 

ys 

76 

77 



dd_ 3 
Sd_ 2 
Sd_i 

<Sc? 2 
8^ 



$d* 



8« 

<5e 2 

& 3 



ce 6 



8/-2 

V-i 

8/o 
8/. 

a/ 3 
8A 

8/e 

8/7 



8y-i 
8y 

8yi 

8y 3 

8?4 

8y 5 
8y 6 

8<ft 



8* 
SA, 

8A 2 
8A 3 
SA 4 

8*6 



So 

^1 £/W 

<$2° 2 2^i 

<5£ 2 ^mj 

0*3 Ct 2 

<5& 3 £m 2 
8»4 S7 8/ 3 
d«4 om 3 

Sis 2i 

* Ms 
ci 6 



And let the differences of V , V,, V 2 , V 3 , &c, be extended in like 
manner and with an analogous notation, only omitting the symbol S. 
Then the complete list of the adjusted fourth differences, expressed 
in terms of the focal quantities fa., <j>i, fa, fa, will be 
d_ i +dd_ i =d_ 4 + 0o 

d_ 3 +M_ 3 =d_ 3 — 80o+ fa 

«?_2 + 8rf_ 2 = ^-2 + 280 o — §01+ 02 

r/_i + &jLi=«Li— 5600 + 280,— 80 2 + 3 

d +Bd =d +7O0 o -560 1 + 280 2 — 80 3 
di +Mi =d } —5600 + 7001—5602 + 2803 
d 2 +id a =dz +2800—5601 + 7002—5603 
d 3 +2d 3 =d 3 - 80o+280i-560 2 + 7O03 



d 4 +Zd t =d t + 0o— 801 + 2802-5603 

d 5 + M& =rf 5 . . . . + 0i— 80 2 +2803 

d 6 +Sd 6 =d 6 + 02— 803 

d 7 +Se?7 =d 7 + 03 

The adjustments under consideration are specifically those of 
the fourth differences d , d u d 2 , d 3 ; but the four supplementary 
differences which precede, and also those which follow, are included 
in this list, in order to avoid any material disturbance of these 
adjoining differences, by operating upon the results which include 
the whole of the effects developed by the focal adjustments. 

Now, in order to determine the focal quantities fa, fa, fa, &c, 
when the sum of the squares of these adjusted differences is a mini- 
mum, the resolving equations, according to the rule on page 145, 
will be found by first multiplying each equation by the coefficient 
of fa, and equating the sum of the products with zero; and then 
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m + $m =0 ' 
m 1 + dm i =0 

m 2 + <S»» 2 = 

m 3 -\-dm 3 =Q 



or 



(») 



severally proceeding to operate in the same manner with respect to 
0i > 02> 03- Hence, observing that 

«L 4 — 8eL 3 + 28rf_ 2 — 56c?_, + 70rf„— 56<A + 28d 2 — 8d 3 +d t 
is the eighth difference of the series 

«-4) <*-3) "-i> **-U <*0» <*1> *2> "3> "4> 

and therefore equal to ?w ; and that the other similar sets of pro- 
ducts reduce in exactly the same manner, the required equations of 
condition are simply 

g»2 =— m ' 
hm,i= — m x 
$m 2 = — m 2 
dm 3 = — m 3 

These show that the twelfth differences of V , V u V 2 , &c, 
must be neutralized by the sixteenth differences of <f> 0) U 2 , &c. 

The sixteenth differences, viz., &w , &Wi, $m 2 , Bm 3) being thus 
known, the values of <j> , <j> u 2 , <j> 3 , are, according to the table of 
coefficients, page 311, vol. xi., determined by the equations 
12870^0—114400!+ 8008fj— 4368<j> 3 =e~m \ 
— 11440^0+ 12870^!— 1144O0 2 + 8OO80 3 =Smi I . ( „. 
8008^0— 1144O0 1 + 1287O^ 2 — 1144O0 3 =gm 2 [ ' ' A ' 
—4368^0+ 8008^— 1144O0 2 + 12870* 3 =gm 3 J 
By solving these equations we get the following formulae for 
computation with four foci : — 

2771340 o = 605&n +1320&M 1 + 1188S»j 2 + 440Sm 3 \ 
2771340 1 = 132Ogw o +3165ami + 3O48gm 2 + 1188^ 
277134f,=1188Sw + 3048Sm, + 3165&W2 + 1320Sm3 
277134^3= 440gOT +1188Sm, + 1320am 2 + 6QbZmJ 

The determination of these coefficients is indeed evidently the 
same as would be obtained by the former method when applied to 
eighth instead of fourth differences ; and they are here similarly 
tabulated up to seven differences, or seven foci. 

Tables B. — Coefficients for deducing Focal Adjustments from Twelfth 

Differences. 



(9) 



12870.0 24310.. 



24310.1 



277134.0 



25 


40 


20 


40 


73 


40 


20 


40 


25 



605 


1320 


1188 


440 


1320 


3165 


3048 


1188 


1188 


3048 


3165 


1320 


440 


1188 


1320 


605 
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33 


88 


108 


72 


22 


88 


253 


328 


228 


72 
108 


108 


328 


449^ 


328 


72 


228 


328 


253 


88 


22 


72 


108 


88 


33 



293930.0 



1859 


5720 


8580 


7800 


4290 


1144 


5720 


18755 


29480 


27780 


15720 


4290 


8580 


29480 


48455 


47480 


27780 


7800 


7800 


27780 


47480 


48455 


29480 


8580 


4290 


15720 


27780 


29480 


18755 


5720 


1144 


4290 


7800 


8580 


5720 


1859 



125970.0 



1183 


4056 


7020 


7800 


5850 


2808 


676 


4056 


14703 


26520 


30420 


23400 


11466 


2808 


7020 


26520 


49695 


58920 


46620 


23400 


5850 


7800 


30420 


58920 


72195 


58920 


30420 


7800 


5850 


23400 


46620 


58920 


49695 


26520 


7020 


2808 


11466 


23400 


30420 


26520 


14703 


4056 


676 


2808 


2850 


7800 


7020 


4056 


1183 



The coefficients of each higher table may, with a considerable 
saving of calculation, be deduced from those of the table imme- 
diately preceding it by the following method : — 

Distinguishing by accentuation the symbols which appertain to 
the preceding table, let the equations analogous to (6) be 

— A^'o+Aq^!— A^'j ±A n _ 2 f n _ 1 =Sm 1 

Aj$' — A^',4- A o 0' 2 zfA n _ 3 <j>' n _i=Sm i 



•+-K-i<l>'o±K-2<i>'i+A H -3t't- • +A f B _l=S»»»- 



(7) 
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where A , Aj, A 2 , &c, are factorial functions of n; and let the 
solution of these equations, as given by the tabular coefficients, be 
denoted thus : — 

<fo =0'o,O^ WJ O +^'<>, 1&»1 +0'o,2^2- • • •+^'o,n-l2»» n _ 1 \ 

0'l =tf>'l,0<K) +<t>'hl% m l +<P'i,-£*»h +t'\.n-\$m n _ 1 I 

$1 =^'2,0^0 +f2,i^ OT i +0'a, aS*w 2 . . . . +0 , 2,«-i^ w *«-i I.. (8) 



0'n-i=f»-i o2»«o + ?>' B _i,iSmi+0' B _,, 2 Sm 2 . . +0'„_i, B _i^m„_ 1 / 

Let the equations of condition which have to determine the 
coefficients of the sought table be 



Ao<Po— Aitf»i + A 2 2 +A n _ 1 ^„_ 1 ±A„^„ =Sm 

— A^o+Ao^! — A^ 2 ±A B _ 2 ^„_ 1 +A n _ ) ^„=gm 1 

A 2 ^. — Ai^! + A ^ 2 : FA B _ 3 ^ n _ 1 ±A B _ 2 0„=g»j 2 

= FA B _ 1 0o±A B _ 2 ^ 1 qrA„_ 3 ^ 2 . . . + A (f> n _ l —A 1 <j>„=dm a _i 



.(9) 



Now, if 8m TA„^ a , 8m,±A n _,0 B , Sm 2 +A„_ 2 n , 

Sm„-i + Aif„ be substituted in place of §m , &»!, Bm 2 , . . . 8»n B _i 
respectively in the equation (7), they will be identical in form with 
the equations (9), omitting for the present the last one. Hence, 
by a like substitution in the first of the formulae (8), we get 

0o=0'o, o(S™o+K<t>«) +$'(,. iO»i±A„_i0 B ) 

+ 9'<>,2(&»2+A„-2f,) +f'o,,-l(H-l+Al^) 

=t'o + (Aip'o, „-i — A 2 0.' Oj B _ 2 +A 3 (p' 0i n _ s —&c.)<j>„. 
And in this way the formula? (8) severally give 



where 



0o =<p'o +^otf'» "» 
<t>i =<f\ +^10„ 
i>2 =^'2 + ^20» 

fn-l=<l>'a-\ + K-l<l>n 



.(10) 



\> = Al^o,n-l ~-A 2 0' OjB _ 2 +A 3 ^' 0jB _ 3 — &c. 
\ t =A,0' ljB -i — A 2 0'i )B _ 2 +A 3 ^'j )B _ 3 — &c. 
X 2 =A 1 ^' 2) „_, — A 2 ^' 2)B _ 2 +A 3 f' %n _ 3 — &c. 

^»- 1 = Ai0'„_ 1, B _ 1 — A 2 0'„_ u B _ 2 + A 3 ^'„_ , f „. 3 — &c. 
= Ai^'o, — A t '<j>' 0i 1 + A 3 0' Oi j— &c. 

The last of the equations (9) now gives 



(") 
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± K(t' o + \>t<d+ K-\{t>\ + Xl0„) ± A„_ 2 (f 2 + X 2 f.) 
.... -A 1 (f„_ I + X„_,^„) + A ^„=aw„. 
But 

±A,.^'oTA„_ 1 ^' I ±A„_ 2 f2 • • • • -Ai0'»_i 

= -A 1 (0'„_ ljO Sw o +0'„_ 1 ,,g>n 1 + 0'„_ Ii2 <5»! 2 +&c.) 
+ A i (<j>'„_ % Sm + f'„_2, iSmi + <j>' n _ 2 , 2 Sm 2 + &c.) 
- A 3 (f„_ 3j gm + 4>'„_ 3j ,3m! + ^>'„-3, 2Sm 2 + &c.) 
&c. &c. 

= — (\ Sm +\ 1 dmi + \ 2 Em 2 .... +X»_i2»»„.i); 
.-. (A -A 1 \„_ 1 +A 2 X„_ 2 . , . . ±A„\o)^>„ 
^X Smo+Xi2mi + X 2 S>W2 .... +X„_iSm„_ 1 + Sm„. 

Hence, as the coefficient of $m„ in the value of ^„ is ^», B =^o,o> 
we must have 

0»=0o,o(Xot>»io+X 1 Sm 1 + X2Sm 2 • • • • +X»_iSm„_ 1 + S»» B ) . . (12) 

where the several coefficients are Ao0o,o> \<f>o,oy ^20o,o> • • • 0o,o> an< i 
these must be the same as the coefficients of <j> , taken in a reverse 
order. But O i g t fle same as the first supplementary difference in 
the scheme, for which a general formula has before been found. 

If p denote the order of the differencing, counted from to the 
differences to be adjusted, then, according to (k), n being now n+p, 

x+1 . . . x+p-1 n+l-x . . . n+p-x . 

<4 ft =2 — — - t*- om. . . (13) 

70 n+p + l...n+2p 2...p-l K ' 

As the coefficients of <j> n are those of <j> reversed, put n—x for 

x, and 

n-x+\ . . . n-x+p-1 x+l...x+p 

<4„=2 £ • — — —cm. .. (14); 

? n+p+l . . . n+2p 2...p-l * K J 

the respective coefficients of which are 

n+1 . . . n+p — 1 



<l>0, n =\>to, 0=P 



_ l =\ l f 00 =p—— 



n+p+l . . . n+2p 

p + 1 «... n+p — 2 
1 n+p + l . . . n+2p 



p + 1 p+2 n—l...n+p — 3 

to, »-2=X2?o, o-P -j Y" n+p+l... n+2p 

p + 1 p+2 p + S n — 2 ... n+p — 4: 
to, „-3=X 3 ^o, o-P -j 2 3- ' n+p + 1 .,. n+2 p 



p p + 1 p + 2 n-2 . . . n+p — S 
to, 3 =K-zto, 0- j- • ~z 3- • n+p + 1 ■ ; ; n + 2p 
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p p + 1 n—l...n+p — 2 



00, 2— A„_20O,O — ' 



2 n+p + 1 ...n + 2p 



p n . . . n +p — 1 
to, F*.-iK.= J * n+^ + l . . . n+2p 



.-.\ = 

«+jt> 

., P + 1 
^i=P 



W+l . . . »+J9 («+j») 2 

n+p + 1 . . . «+2p~~ n(n+2p) 
P 



i.o= — ■ , , — •— ™r^ — rzrroi^'o.o . . . (15) 



X 2=P 
X 3 =p 



1 (n+jt>)(»+p-l) 
p+1 p + 2 w(w-l) 



1 2 (n+p)(n+p-l)(n+p-2) 

p+l p + 2 p + 3 n(n-l)(n-2) 



\i-3== 



12 3 (»+jB)(»+p-l)(»+/>-2)(w+/>-3) 

jD p+1 p+2 n(n-l)(»-2) 

T ' ~~2 3~~ ' (» +i>)(w +p - 1 )(n +p - 2) 



_jo p+1 n(»-l) 

A «-2 — T " " 



1 2 («+p)( w +p-l) 
/> *» 

A„_i=r • — — 

1 n+p 

K =1 (16) 

These values may be conveniently made to depend successively 
upon each other, as follows : — 

K =1 X„=-^-X„ 

n+p 

x - P+ l > 

P + 2 



A.. 


-i— 


T A 




a„. 


•2= 


»-l 

"2" 


X. 


A„. 


-3= 


«-2 
"IT 


*2 



«+p-2 
_ Z +3_ 

&c. &c (17) 

The coefficients are hence all easily deduced by means of the 
formulae (12) and (10). 

Moreover, according to the formula (10) the absolute value of 
the function $ x is obviously S(A a .^ M ), in which the terms under 2 
are determined by giving to n the descending series of values n, 
n— 1, ra— 2, &c. Hence from the foregoing general expressions 
for 0„ and the factors A, we are thus enabled to obtain, for the 
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independent and direct calculation of any focal coefficient, the 
following general formula: — 

n — x+1 . . . n—x+p n—y+1 . . . n—y+p 
*'■ "~ 2...p-l 2777^1 X 

„ (x+l .. .x+p-l)(y+l .. .y+p-1) ^ _ ^ 
n + 1 . . . n + 2jb 
where the terms under S are to be found by successively diminish- 
ing by unity every factor in both numerator and denominator. 
The number of these terms after the first will evidently be the 
lesser of the two values x, y. Also if x'=n— x, y'=n—y, be com- 
plementary numbers of the table, since then <j> Xty =<j> x ' ty ', the calcu- 
lation is simplified, and consists of less terms, if we choose from 
these two equivalent forms that which contains the least of the four 
values x, y, x', y'. 

The formula (18) expresses the value of the focal coefficient 
which belongs to a table constructed for n+1 differences. Let 
this, for the present, be distinguished by placing the symbol n as 
an ordinal exponent, and not as a power; thus <j>^ y designates the 
y + lth coefficient of the x + lth line of a table for n + 1 values. 
Then, by applying the formula to express the several values of the 
following coefficients, viz., 

*:„. Ct„-.. C:l-2> — c:;,o> o, 

the leading factors, which precede 2, will remain unaltered, but each 
succeeding series under 2 will become divested of its first term. 
The successive terms of the formula, taken separately, must, there- 
fore, be the respective values of the first differences or decrements 
of the stated series of functions. Hence, by dividing the first of 
these terms by the second, we find that 



:) 



••• ♦;,=^L,-i+ jqp^ — - x -y — (*.-!.,-.-*.-*,-»)• -09) 

which supplies an easy method of calculating a coefficient for any table 
from those contained in the two tables immediately preceding it. 
In extreme and exceptional cases, we have 



/. n x+p — 1 p\ . , , 



n-\-2p i 

n-\-p x-\-p — 1 



x-1,0 



A »-l 



n+2p x r *- , -° 



_n+P_ n+p 
9o -° n+2p' n 



(20) 
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For tables A, p=4* ; for tables B, p=S. 

Examples. — In the table B for seven values, we have 

, . , 6 11 9 , , . . 
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22 4 2 KY3 - 1 
2778 27 / 2778 



29393 

10878 



+ 



8 V29393 
1554 



54 \ 
4199/ 



46620 Y 
12597oJ' 



29393 4199 \ 
. ( ■. . 6 11 8\ , 

780 / _ 23400 "\ 



=7 



780 



29393" 4199N. 125970/' 



A 8 



14 11 

= 22* 4 " 

7 780 



4 29393 



195 / 5850 \ 
: 4.199 V~ 125970./' 



Also in relation to the tables B, the computed values of <p 0tQ 
A , Xi , &c., for a few values of n are here appended : — 



n 


^0,0 


*o 


Xi 


x 2 


X3 


K 


\5 


K 


X 7 


*8 


x 3 


X10 


1 


9 


8 
9 


1 




















17.1430 


2 


5 


4 

5 


8 
5 


1 


















17.286 


3 


55 


8 


108 


24 


1 
















17.19.78 


11 


55 


11 
















4 


33 


2 


24 


36 


8 


1 














17.19.26 


3 


11 


11 


3 
















143 


8 


30 


600 


60 


40 














5 


17.19.70 


13 


13 


143 


13 


13 


1 












6 


91 


4 


216 


450 


600 


540 


24 


1 










17.19.30 


7 


91 


91 


91 


91 


7 












195 


8 


12 


72 


110 


120 


36 


56 










'/ 


17.19.46 


15 


5 


13 


13 


13 


5 


15 


1 








8 


130 
17.19.23 


1 

2 


12 
5 


6 


132 
13 


165 
13 


12 


42 
5 


4 


1 








17.26 


8 


81 


108 


198 


12.297 


297 


252 


162 


72 






9 


19.23.45 


17 


34 


17 


17 


13.17 


17 


17 


17 


7 


1 




10 


17.18 


4 


40 


225 


220 


330 


396 


385 


300 


180 


40 




19.23.25 


9 


17 


24 


17 


17 


17 


17 


17 


17 


9 


1 



N 2 
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We shall conclude with a summary recapitulation of the general 
principles of calculation, according to the method propounded in 
the last investigation. 

Let the given series of numbers, about the vicinity of any 
objectionable irregularities amongst the fourth differences, be differ- 
enced partially out to the twelfth order of differences. For example, 
in the following scheme, as in that on page 164, four refractory 
fourth differences (d , d u r/ 2 , d 3 ) are considered to be centres of a 
disturbance which demands adjustment. 

«-3 



v, 
v 2 
v 3 
v 4 

v 6 
v 7 



a 
a 2 

«3 

o 4 

«7 



"-1 

»1 

b 2 

h 
h 

b, 



C-3 

C-2 

e_] 

Co 

C\ 
r 2 
e 3 

°1— 
C 5 

C6 

c 7 



d-3 
«L 2 

A 1 

d 
d x 
d 2 

A. 

d t 

d 5 
de 
d-j 



/-. 

/-l 

e -l 

— /o 

*0 

/. 
/. 
/. 

A 

■A_ 
A 
A 



9-1 


h a 


yo 


h, 


9i 


A 2 


9-2 


A3 


i/S 


K 


9* 


h 


9b 


h 


96 


hf 


97 





h 

*2 
*3 
K 

h 



m 
m 2 



Thus, by carrying out the differencing, the four twelfth differences 
m , mi, m 2 , m 3 are found ; and after reversing their algebraic signs, 
according to equations (w), they will then become the values of 
8m , Swij, §m 2 , Swig. 

Now, by table B for four differences, the focal values are 

605S>w + 1 320Smi + 1 1 88Sm 2 + 44Qgm 3 
^ 0— 277134 

_ 1320 3m +3165gOTi + 3048g?w 2 +11882m 3 
^ 1_ 277134 



1188gmo + 3048gOT 1 + 3165gTO 2 -H320gOT 3 
277134 



440Sm + 1 188am! + 1320Sm 2 + 605fo« 3 
277134 



Having calculated these numerically, their expanded fourth differ- 
ences are then to be developed as in the scheme on page 164, and 
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thus the sought values (y , . . . v 7 ) of the corrections or adjustments 
of eight primitive numbers (V , . . . V 7 ) are determined. 

These adjustments will affect the values of twelve fourth differ- 
ences (rf_ 4 , . . . a\), which thereby become so modified and reduced 
in magnitude that the sum of their squares is less than would have 
similarly resulted from any other values of the focal quantities (0 O > 
4>i > i>2> <p3)> an d ft t the same time the values of the principal 
adjustments (v , . . . v-,) are such that the sum of their squares is 
less than that of any other set of numbers that can produce, as 
fourth differences, the four middle adjustments (§d , Bd u 8rf 2 , Sd 3 ). 

If, according to the same method, second differences be used, 
instead of fourth, the operations will be perfectly analogous, but 
will, of course, involve less calculation. In the corresponding case 
of four centres, the schemes will be as follows : — 



V n 


«U 


v, 


«0 




a, 


v, 




v 3 


a 2 




a 3 


M 






a 4 


»s 






«5 



i-2 




*-, 


O-l 


*o 


c 




c l 


i>, 




h 


H 


h 


H 




c 4 


*4 






«s 


h 





/o 

/ 

A 

fz 



¥o=~fo 
2/2= -ft 









85 o 










" 1 




8« , 






v a 




M_, 





H> 


w 


*>i 


a 


8A„ 


8c 


8fl?, 


IV, 




a l 




^! 






»? 




8A, 




&/, 


tv 2 




a 2 




8c 2 






»8 




**l 




Wj, 


te, 




a 3 




8c 3 






»4 




**» 




H 


w A 




«4 




8c 4 






»« 


a 5 




Sc 5 


H 



8e 
&, 

8e 2 

8e 3 
8« 4 



»/o 
V. 

8/3 



and the formula? for finding the focal values from the reversed sixth 
differences will be, according to the table A for four values, viz.: — 



¥>o 



_ 498/ + 848/ + 708/, + 288/ 3 

462 
= 78/ +128/ + 108/ 2 +48/ 3 

6G 
_ 848/ + 1 778/ + 1 648/, + 708 / 3 
462 
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7Qg/ + 164S/. + 177S/ 2 + 84g/ 3 
* 2= 462 ' 

28S/ + 708/, + 848/ + 49S/ 3 
* 3= 462 

_ 48/ +10g/, + 12a/ 2 +7g/3 

66 

Whether any real advantage is derivable from this particular 
method at all commensurate with the labour and complication of the 
process, or whether the method itself admits of material improve- 
ment by the substitution of a different hypothesis, must for the 
present be left to be determined hereafter. 

The former method, or that of effecting, under the best possible 
conditions, any prescribed small adjustments of fourth differences, 
is much more simple in its application, the focal values being 
found at once by using the coefficients contained in the first set of 
tables. 

We may also here state, that the two practical rules given at 
the commencement, for the adjustment of absolute or prominent 
errors, are at all times useful for the verification and revision of 
consecutive calculations, and are indispensable for the preliminary 
adjustment of a series of numbers previous to interpolation, or as a 
preparatory for further adjustment by one of the methods just 
alluded to, if this be considered necessary. 

Should the resulting adjustments happen to be of any con- 
siderable magnitude, the fact will indicate either that the assumed 
interval of the table is too great, or that the presumed approximation 
of the given series of numbers to the true values is not entitled to 
be much relied upon. 

In verifying results of computation, it should be borne in mind 
that, though practically correct, certain slight deviations from the 
exact values are continually incurred by the necessity of cutting off 
the figures at a given place of decimals. The last digit being 
always put down at the nearest figure, the amount of inaccuracy 
from this cause can never exceed half an unit of the last digit, and, 
amongst the several values, it may be either in excess or defect. 
Hence, if three consecutive differences, according to the rule, indi- 
cate an error of only an unit in the last figure, the case may in 
general be passed over without suspicion of inadvertence, since a 
displacement to that extent may naturally be accounted for by the 
occasional occurrence of something like half an unit alternately in 
excess and defect of two consecutive numbers; but if the discrepancy 
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of the differences should indicate a primary displacement of more 
than one unit, we may then be assured that an error really exists. 
The practical observance of this rule will obviate any unnecessary 
and useless re-examination of the original calculations. 

It may be right also to observe generally, that the investigations 
that have been given in this paper have an implied reference to the 
elimination of casual errors or incidental imperfections, and do not 
necessarily apply to errors, assignable to certain specific causes, 
which may admit of being expressed by some mathematical function; 
for in such case the variations involved are no longer promiscuous 
in their character, and they do not contribute to produce any 
marked irregularities in the progression of the tabular differences. 
Such, for instance, are continued errors of principle in calculation, 
and, in some cases, personal errors of observation, which undoubt- 
edly should be obviated or detected independently by other means, 
without having recourse to any theory of differences. 

Note. — After finding the supplementary fourth difference adjustments by 
the formulas (k), an analytical expression for the corrected functions may be 
obtained by completing the scheme of differences of the original functions 
(V , ^li ■ • • V»)> tnat i s > by treating the previous functions as zeros, and 
thus prefixing the values of the corresponding supplementary differences 
(d. t , d_ 3 , d_2, d_i), which are simply deduced by numerical subtraction. 
For then, by the formulae (e), the value of an adjusted function will be 

v m+ * m=S ' (»-^X»-y»X»-i+ i) R+R)> 
or { y +v)m ^ m -^ 3 ^-f'^ m -^ %i + M%, 

in which the ordinal characteristic (x) begins with the antecedent value —4. 
Consider the case in which the fourth differences are absolutely neu- 
tralized, or in which dd = — d , Srfj= — d lt &c. Then, in developing the 
terms, under 2', of this expression, all of them, after the first four, will 
vanish; because the adjusted fourth differences d t -\-M x become severally 
zero when x reaches the values 0, 1, 2, &c. Therefore, in this case, 

+ (»-lff»+l ) (rf+M) _ i 

(m— 2)(m — Vim,, „ ,. 
+ - £3 J ~( d +MU 0); 

and this algebraic function of the third degree in m is that which represents 
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the original series of values (V , V^ . . . V„) more approximately than any 
other algebraic function of the same degree. 

And the most approximate algebraic function of any degree may be 
obtained in like manner, by employing differences of an order relatively one 
degree higher. The analogous geometrical problem -would be the deter- 
mination of that curve line of a given order which shall follow most closely 
upon the devious track indicated by a number of given points. 



CORRESPONDENCE. 



THINGS WORTH NOTING. 

To the Editor of the Assurance Magazine. 

Sir, — I regret to have observed of late a great falling off in that which 
I consider by no means the least important or least interesting department 
of the Magazine — I mean the Correspondence department. I feel sure 
that if those of us who are accustomed to derive information and instruction 
from your pages were duly to discharge the duty we owe to the Magazine 
in return, the department to which I refer would speedily show signs of 
renewed vitality. We all, in the course of our reading, meet from time 
to time with " Things worth noting," which although often apparently trivial 
in themselves, are yet frequently by no means destitute of interest, either 
from their suggestiveness, their throwing light on some point in the history 
of our science, their illustrating the peculiarities and comparative advan- 
tages of different methods of investigation, or the like. Now, why should 
we not offer, now and then, for editorial approval, if we have nothing 
better, a little "Budget" of such "Things" of the above or cognate character 
as have come under our notice? It is very conceivable that wholesome 
discussion might thereby not unfrequently be excited, and no small amount 
of instruction elicited. 

The advantage to be derived by the adoption of the course I have ven- 
tured to suggest, would not, as all experience proves, be entirely, or even 
principally, on the side of the readers of the Magazine. The writers would 
largely participate. This is a trite point. Suffice it to say, that the infor- 
mation we acquire as to the amount and the completeness of our acquaintance 
with a subject when we try to write upon it, and the stimulus given to our 
endeavours to supply the deficiencies which ^fhen almost always come to 
light, furnish far more than a compensation for any amount of labour that 
the effort to put our ideas upon paper may have cost us. 

" Example," we are often told, " is better than precept." I thereforo 
propose here to exemplify, by a first instalment, what I mean by " Things 
worth noting." 

1. The history of the celebrated formula, a x =vp s ,(l+a ll+i ), which 
assigns the value of an annuity upon (x) in terms of that of the corresponding 
benefit upon (x+ 1) is pretty well known, up to a certain point. Mr. Milue 
("Introduction," pp. xv., xvi.) attributes its discovery to Thomas Simpson, 
who published it in 1742, and he further informs us that it was rediscovered 
by Euler, whose publication of it dates in 1760. Subsequent writers, as 



